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Abstract: Let g > 2 be a positive integer, i? be a fractional Brownian motion with Hurst index 
H G (0,1), Z be an Hermite random variable of index q, and Hq denote the Hermite polynomial 
having degree q. For any n > 1, set Vn — J2'k=o Hqi^k+i — B^)- The aim of the current paper is 
. to derive, in the case when the Hurst index verifies H > 1 — l/(2q), an upper bound for the total 

I variation distance between the laws ^(Z„) and ^{Z), where Z„ stands for the correct renormal- 

^SJ ' ization of Vn which converges in distribution towards Z. Our results should be compared with those 

\—i . obtained recently by Nourdin and Peccati (2007) in the case when H < 1 — l/{2q), corresponding to 

the situation where one has normal approximation. 
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^ ; 1 Introduction 

Let g > 2 be a positive integer and i? be a fractional Brownian motion (fBm) with Hurst 
index H € (0, 1). The asymptotic behavior of the g- Hermite power variations of B with 
^ , respect to N, defined as 

(N. Vn = y^Hg{Bk+i-Bk), n>l, (1.1) 

I has recently received a lot of attention, see e.g. [9], [10] and references therein. Here, Hq 

' stands for the Hermite polynomial with degree q, given by Hq{x) = {—lYe^ '^^^^(^"^ '^^)- 

^0 , We have H2{x) = — 1, H^{x) = x^ — Sx, and so on. The analysis of the asymptotic 

behavior of (jl.ip is motivated, for instance, by the traditional applications of quadratic 
variations to parameter estimation problems (see e.g. [HIHEKIS] and references therein). 
^ ' In the particular case of the standard Brownian motion (that is when if = ^), the asymp- 

totic behavior of (jl.ip can be immediately deduced from the classical central limit theorem. 
When H ^, the increments of B are not independent anymore and the asymptotic behav- 
ior of (jl.ip is consequently more difficult to reach. However, thanks to the seminal works of 
Breuer and Major |3], Dobrushin and Major [6], Giraitis and Surgailis [7| and Taqqu |12j . 
it is well-known that we have, as n — > oo: 
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1. If < if < 1 - l/{2q) then 



2. If if = 1 - l/{2q) then 



^n:=— ^ ^(0,1). (1.2) 



:= ^ ^(0, 1). (1.3) 

CTq^Hyn logn 



3. If if > 1 - l/(2g) then 



Zn '■= — -n — F7T Z ~ "Hermite random variable". (1-4) 



Here, Cij^ff > denotes an (explicit) constant depending only on q and H. Moreover, the 
Hermite random variable Z appearing in (jl.4p is defined as the value at time 1 of the Hermite 
process, i.e. 

Z = I^{L,), (1.5) 

where 1^ denotes the g-multiple stochastic integral with respect to a Wiener process W , 
while Li is the symmetric kernel defined as 

1 

Li{yi,. . . ,yq) = —\Q^i\q{yi,. . . ,yq) I diKuiu^yi) . . .diKH{u,yq)du, 

with Kh the square integrable kernel given by (|2.2|) . We refer to [TU] for a complete discus- 
sion of this subject. 

The exact expression of the distribution function of Z^ is very complicated when H ^ ^. 
That is why, when n is assumed to be large, it is common to use ()1.2p ~ (ll.4p as a justification 
to replace, in any computation involving it, the distribution function of Zn by that of the 
corresponding limit. Of course, if one applies this strategy without any care (in particular, if 
one has no idea, even imprecise, of an error bounds in terms of n), then it is easy to imagine 
that the obtained result could be very far from the reality (as a "concrete" example, see (II. 6p 
below, which represents the worst case we will obtain here). To the best of our knowledge, in 
all the works using p.2p - p.4p with statistical applications in mind (for instance [H HI El 13] ) j 
never their author (s) considered this problem. The current paper, together with [TT], seem 
to be the first attempt in such direction. 

Recall that the total variation distance between the laws of two real-valued random 
variables Y and X is defined as 

dTv {^{Y),^{X)) = sup I P{Y £ A) - P{X e A)\ 

Ae.«(R) 

where ^{M) denotes the set of Borelian of R. In [llj, by combining Stein's method with 
Malliavin calculus (see also Theorem 11.31 below), the following result is shown: 
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Theorem 1.1. If H < \ — l/{2q) then, for some constant Cq^n > depending uniquely on 
q and H , we have: 



n 



tfH€{0,^] 



2' 2i}-2J 

2g-3 . 

2g-2'^ 2q^ 



for Zn defined by 

Here, we deal with the remaining cases, that is when H ^ [1 — ^,1). Our main resuh is 
as follows: 

Theorem 1.2. 1. If H = \ — \/{2q) then, for some constant Cq^n > depending uniquely 
on q and H, we have 

dTv{^{Zn),^iO,l)) < (1.6) 

for Zn defined by il.3]) . 

2. If H € (1 — 1/(2(7), 1) then, for some constant Cq^n > depending uniquely on q and 
H , we have 

drv {^{Zn),^{Z)) < Cq,H n'^^-" (1.7) 
for Zn and Z defined by ( |j.^| ). 

Actually, the case when H = 1 — 1/ {2q) can be tackled by mimicking the proof of Theorem 
ll.li Only minor changes are required: we will also conclude thanks to the following general 
result by Nourdin and Peccati. 

Theorem 1.3. (cf. Fix an integer q > 2 and let {/n}n>i be a sequence of S^®'^ . Then 

we have 

drv {^{Iq{fn)),-yi'i0^)) <2^ E (^l-^\\DIq{fn)\\l^ , 

where D stands for the Malliavin derivative with respect to X. 

Here, and for the rest of the paper, X denotes a centered Gaussian isonormal process on 
a real separable Hilbert space ^ and, as usual, ij®^ (resp. Iq) stands for the gth symmetric 
tensor product of 9) (resp. the multiple Wiener-Ito integral of order q with respect to X). 
See Section [2] for more precise definitions and properties. 

When H ^ {\ — l/(2g),l), Theorem 11.31 can not be used (the limit in p.4p being not 
Gaussian), and another argument is required. Our new idea is as follows. First, using the 
scaling property (|2.1|) of fBm, we construct, for every fixed n, a copy Sn of Zn that converges 
in L^. Then, we use the following result by Davydov and Martynova. 

Theorem 1.4. (cf. J^; see also J^) Fix an integer q > 2 and let f E \ {0}. Then, for 
any sequence {fn}n>i C ^®'^ converging to f, their exists a constant Cgj, depending only 
on q and f, such that: 



drv (^(I,(/„)),^(I,(/))) < CqjWfn - f\\]il. 
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The rest of the paper is organized as follows. In Section [21 some preliminary results on 
fractional Brownian motion and Malliavin calculus are presented. Section [3] deals with the 
case iJ G (1 — l/(2g), 1), while the critical case H = 1 — l/{2q) is considered in Section [H 

2 Preliminaries 

Let B = {Bt, i > 0} be a fBm with Hurst index H G (0, 1), that is a centered Gaussian 
process, started from zero and with covariance function E{BsBt) = R{s,t), where 

R(s, t) = ^ {t^^ + s^^ -\t- s\^^) ; s,t>0. 

In particular, it is immediately shown that B has stationary increments and is selfsimilar of 
index H. Precisely, for any h,c > 0, we have 

{Bt^h-Bh,t>0}''^ {Bt,t>0} and {c"^ S^, t > 0} {S^, t > 0}. (2.1) 

For any choice of the Hurst parameter H £ (0,1), the Gaussian space generated by B 
can be identified with an isonormal Gaussian process of the type B = {B{h) : h G Sj}, where 
the real and separable Hilbert space ^ is defined as follows: (i) denote by the set of all 
M- valued step functions on [0, oo), (ii) define ^ as the Hilbert space obtained by closing 
with respect to the scalar product 

(l[0,t],l[0,s])j5 = R{t,s). 

In particular, with such a notation, one has that Bt = i?(l[Q_^]). 

From now, assume on one hand that B is defined on [0, 1] and on the other hand that 
the Hurst index verifies H > ^. The covariance kernel R can be written as 

rsAt 

R{t,s)= / KH{t,r)KHis,r)dr, 
Jo 

where Kh is the square integrable kernel given by 

KH{t,s)=r(^H+^-^ \t-s)''-'^F(^H-l,^-H,H + l,l-i}j, (2.2) 

F{a, b, c, z) being the Gauss hypergeometric function. Consider the linear operator from 
<r to L2([0,1]) defined by 

{K*H^){s) = Kh{1,s)^{s) + r (v9(r) - ^{s))diKH{r,s)dr. 

J s 

For any pair of step functions ip and tp in S, we have {K'^cp, K'^ijj) = (99, -0) $5- As 
a consequence, the operator K'^ provides an isometry between the Hilbert spaces and 
L^([0, 1]). Hence, the process W = (Wt)tG[o,i] defined by 

Wt = B{{K*HrHl[o,t])) (2-3) 
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is a Wiener process, and the process B has an integral representation of the form Bt = 
KH{t,s)dWs, because l[o,t])(s) = KH{t,s). 

The elements of ^ may be not functions but distributions. However, S) contains the 
subset of all measurable functions / : [0, 1] ^ M such that 

\f{u)\\f{v)\\u - v\^^-^dudv < oo. 

[0,1]2 

Moreover, for f,g£ \S^\, we have 

(/, g)s^ = H{2H - 1) / f{u) g{v) \u - v\''"-^dudv. 

i[0,l]2 

In the sequel, we note and i^®'', respectively, the tensor space and the symmetric 
tensor space of order q>l. Let {e^ : /c > 1} be a complete orthogonal system in i^. Given 
/ G and g G ^5®"^, for every r = 0, . . . ,p A q, the rth contraction of / and g is the 
element of i^QCp+'J^^r) (jgfined as 



f®rg= ^ (/,eji (g) • • • (g)ei^)fl®r (g) (5f,eii 



In particular, note that f ®q g = f ® g and, when p = q, that f ®p g = {f,g)^. Since, in 
general, the contraction f ®rg is not a symmetric element of ^'^(p+'?-2r) ^ define f®rg as 
the canonical symmetrization of / ®r g- When / G •^®'', we write Iq{f) to indicate its gth 
multiple integral with respect to B. The following formula is useful to compute the product 
of such integrals: if / G ij®'' and g G ^®'^, then 

Whig) = ) (^) W2r(/ir<7)- (2.4) 

Let =5^ be the set of cylindrical functionals F of the form 

F = ^{B{hi),...,B{hn)), (2.5) 

where n > 1, hi G and the function ip G '^°°(]R") is such that its partial derivatives have 
polynomial growth. The Malliavin derivative DF of a functional F of the form p.5p is the 
square integrable i^-valued random variable defined as 

n 

DF = ^dMB{hi), . . . ,B{hn))h,, 
1=1 

where dup denotes the ith partial derivative of ip. In particular, one has that D^Bt = l[o,j] {s) 
for every s,t G [0,1]. As usual, D"*^'^ denotes the closure of =5^ with respect to the norm 
II ■ 111, 2; defined by the relation ||-F||f 2 — ^\-^\ + -^11-^ -^11^- Note that every multiple 
integral belongs to D^'^. Moreover, we have 

A(/,(/))=gVi(/(-,t)) 
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for any / € S)®'^ and t > 0. The Malliavin derivative D also satisfies the following 
chain rule formula: if ip : M" — > R is continuously differentiable with bounded deriva- 
tives and if {Fi, . . . , Fn) is a random vector such that each component belongs to D^'-^, then 
^{Fi, . . . , Fn) is itself an element of D^'^, and moreover 

n 

D^{Fi, ...,Fn) = Y^ dMFi,. . . , Fn)DFi. 



3 Case H e {1 - l/(2g),l) 

In this section, we fix (7 > 2, we assume that H > 1 — and we consider Z defined by (II. 5p 
for W the Wiener process defined by (12. 3p . By the scaling property ()2.ip of fBm, remark 
first that Zn, defined by (|1.4p . has the same law, for any fixed n > 1, as 

n-l 

Sn = n^(^-^)-i - B,/n)) = Wn), (3.1) 

fe=0 

for fn = n'^^^ X]fc=o '^fk/n (fc+i)/n] ^ LlflJ) Theorem 1 (point 3), it is shown that the 

sequence {5„}„>i converges in towards Z, or equivalently that {/n}n>i is Cauchy in 
ij®"^. Here, we precise the rate of this convergence: 

Proposition 3.1. Let f denote the limit of the Cauchy sequence {/n}n>i in S^®'^ . We have 
E\Sn - Z^ = E\lgifn) - /,(/)!' = Wfn " f\\le, = Oin^'^-''^'^'' ) , as 00. 

Proposition Em together with Theorem II. 41 above, immediately entails p.7p so that the rest 
of this section is devoted to the proof of the proposition. 

Proof of Proposition lg.il We have 

n-l 

ll/nlll©, = "-^"^"^ ^ (l[fc/n,(fe+l)/n], l[Z/n,(«+l)/n])|, 
k,l=0 

n-l / i-{k+l)/n r{l+l)/n 

kl=0 yk/n Jl/n J 



By letting n goes to infinity, we obtain 

'2 - Hi{2H-l)'' [ \u-v\'^'i^-^'^dudv 

J [0,1? 



n-l 



"--^ Mk+l)/n /.(/+l)/n 

= H^{2H-iyY, / / dv\u-v\'^'i^-^^. (3.3) 

Now, let (j) £ 1^1 . We have 



n-l 



1=0 
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n-l / .a+l)/n /■! ^ 
' 2H-2 



1=0 -^0 



dv / du 4>{u)\u — v\ 



By letting n goes to infinity, we obtain 

9 



\2H-2 



Hence, we have 

/.I / r(k+l)/n 

fc=0 



k^oJo \Jk/n 

r(l+l)/n ( f 

= H''{2H-lfn''-^ / 



k/r, 

( /■(fc+l)/n 

du\u-v\'''-'' 1 . (3.4) 



'■(«+l)/n / /-(fc+ll/n \ 

2H-2 



kjn 

Finally, by combining (j3.2p . (jS.Sp and (|3.4p . and by using among others elementary change 
of variables, we can write: 

f / i-(k+\)ln r{l+l)/n \'i 

\\fn-f\\%o. = H^ilH-lf Y^ln"^-^ {j^ duj^ dv\u-v\^''~^\ 



k,l=0 



r.(«+l)/n / f(k+l)/n ^ 

-2n5-i / dv{ du\u-v\^^-^ 



l/n \Jk/n 
fik+l)/n f{l+l)/n 

+ dw dz\w - zl'^i^-^'^ 

J k/n J l/n 



H'>{2H-iyn^i-^-^i'' V If / du [ 



dv\k-l + u-v\'^"''^ 



k,l=0 

-2 I dv(^j^ du\k-l + u-v\^^-'^^ + j^duj^ dv\k-l + u-v\'^'^"-^'^ 



\Jo Jo 



1 rl \1 
dv\r + u- 



re 

1 rl 

2qH-2q 



-2 I dv( I du\r + u-v\'^^ "A + I du I dv\r + u - v 
Jo \Jo J Jo Jo 

Consequently, to achieve the proof of Theorem 13. H it remains to ensure that the sum over 
Z in (13. 5p is finite. For r > 1, elementary computations give 



(3.5) 



(^^ duj^ dv\r + u-v\^''-^^ = (2i7(2F-l))-^((r + l)2^-2r2^ + (r-l)2^)'' 

= (^^2H-2^Q^^2H-A^y 
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and 



1 .1 
du / dv\r + u — v\ 

^0 



2qH-2q 



_^ -|^^2gH-2g+2 _ 2^2gi?-2g+2 _^ _ i-j2qH-2q+2 



{2qH -2q + l){2qH - 2q + 2) 

^2qH-2q _^ Q^^2qH-2q-2^ _ 



(3.7) 



Moreover, using the inequality |(1 + x)^^ ^ — 1 - {2H — l)x\ < {2H — 1){H — for 
X S [0, 1], we can write 



dv I / du\r + u — V 



2^-2^ = (2H - 1)-' l\{r + 1 - vf""-' 



[r-vf'-ydv 



i2H - j\r - vf^^~^ ((1 + ^)^^-^ - l)' 

I {r-vf 
Jo 



dv 



r — V r — V 



where the remainder term R verifies |-R(u)| < (1 — H)v?' . In particular, for any v € [0, 1], we 
have 



(r — v) 



Ri 



1 



"r — V 



< 



1-H 
r — 1 



Hence, we deduce: 



dv ( / du\r + u — V 



2H-2 



{r - vfi^-'^'^ {1 + 0{l/r))Uv 

1 _ n _ 1 /r\'2qH-2q+l 

= ^2gH-2,+i l (1 yr) (l+0(l/r)) 

2qH-2q + l ^ ^ ' " 

= r^ili-^i+\i/r + 0(l/r2)(l + 0(l/r)) 

^ ^2qH-2q^(j(^^2qH-2q-ly (33) 

By combining (j3.6|) . (|3.7|) and (j3.8|) . we obtain (since similar arguments also apply for 
r < -1) that 



1 ^1 \ 9 

du / dv\r + u — vl"^^ '"^ 
Jo 



dv ( / du\r + u — V 
-1 /■! 



2H-2 







+ / du dv\r + u — v\ 
Jo Jo 



2qH-2q 



is 0(|rp'^^ SO that the sum over Z in ()3.5p is finite. The proof of Proposition 13.11 is 

done. □ 



4 Case H = 1-1/ {2q) 

As we already pointed out in the Introduction, the proof of (|1.6|) is a slight adaptation 
of that of Theorem 11.11 (that is Theorem 4. 1 in [IT] ) which was devoted to the case when 
H < 1 — l/{2q). That is why we will only focus, here, on the differences between the cases 
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H < 1 — l/{2q) and H = 1 — \/[2q). In particular, we will freely refer to [TT] each time we 
need an estimate already computed therein. 

From now, fix if = 1 — \/{2q) and let us evaluate the right-hand side in Theorem 11.31 
Once again, instead of Z„, we will rather use Sn defined by 

in the sequel, in order to facilitate the connection with [TT]. First, observe that the covariance 
function pu of the Gaussian sequence — B^j^fj^yQ, given by 



pH[r) = -{\r + 1|2-V<? _ 2|r|2-V. + |^ _ 



verifies the following straightforward expansion: 

PHirf = ((1 - ^)(1 - ^)) Vr' + O(kr'), as \r\ ^ oo. 
Using n2"^/^(l[fc/„,(fc+i)/„], l[i/n,(i+i)/n]>ij = PH{k - I), note that 



(4.1) 



Var(5„) - ^2_,^^„ X] -^[■^5(-'-^/n,(fc+l)/n])^'?(^^/n,{«+l)/n])] 



n— 1 I n— 1 



^ k,l=0 H b fc;=o 

from which, together with (j4.ip . we deduce the exact value of a^^: 

:= hm V - If = 2ql ( ' . (4.2) 

In order to apply Theorem II. 3^ we compute the Malliavin derivative of Sn'- 

_ qn s-^ \ -, 



Hence 



0-2 loffn ^ ^'l~^i'^tk%\k+l)/n])^'l-'^i^tlJri}l+l)/n])(^[k/n,{k+l)/n]'^^ 
H ^ kl=0 



The multiplication formula ()2.4p yields 



X 
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1 

r+1 

k,l=0 



X Yl ^2g-2-2r(l^/„^^(fc+i)/„]®l^/„^(Vl)/n])(l[fe/n,(A:+l)/n]i 1 

We can rewrite 



l-h\DSn\\l = l-YAr{ 



q-1 

' . [n] 

r=0 

where 

2 

, / a- i ^ 
rr! 



ajj log n 

n-l 

X ^ ^2g-2-2r(l§/„^(fc+i)/„]^lf^/\(;/l)/„])(l[fc/n,(fc+l)/n], l[Vn,(«+l)H 
fc,Z=0 

For the term Ag_i(n), we have: 

n-l 



l-Aq_i{n) = I - ^^^^ ^^^ ^ ^ (l[fc/«,(fc+l)/n],l[;/n,a+l)/n])^ 

/c,/=0 



^2 iQgj^'^ / > \-L'>./'i.A'>-T~J-;/"I' ~L'/"A'-i-J-;/"J'f) 



n-l 



^ p^(A; - = 1 - Y.{n- \r\)pH{rY 

-'^ fc,/=0 -'^ ^ |r|<n 



1 - E '«^('')' + E IH/'^^W' = 0(1/ log n) 

'"^^1°^^!^^^ a^nlogn^^^^ 



where the last estimate comes from the development (|4.ip of pn and from the exact value 
(lOI) of al^. 

Next, we show that for any &xed r < g — 2, we have E\Ar{n)\'^ = 0(1/ log n). Indeed: 

E\Ar{n)\'' 



- c{H,r,q)- — ^ (l[fc/n,(A:+l)/n], l[«/n,(/+l)/n])^'*'^(l[fc/n,(fe+l)/n], l[//n,{«+l)/n])^"''^ 

i,j,k,l=0 

^ \-^[fc/n,(fc+l)/n]'<2>'-^[«/n,(/+l)/n]' ■^[j/n,(i+l)/n] J-[j/„,(j+l)/n] /i3®29-2-2r 

= X] X] c{H,r,q,a,P,-f,d)Br,a,p,'y,5{n) 

7,<5>0 a,/3>0 
7+(5=g— r— 1 Q!+/3=g— r— 1 

where c(-) is a generic constant depending only on its arguments and 



4g-4 



n-l 



-Br,a,/3,7,5(n) = -j k 

logn^ 

i,j,k,l=0 



(l[fe/n,(fe+l)/n], l[//n,{«+l)/n])^^^(l[Vn,(i+l)/n], l[i/n,(j+l)/n])^'*'^(l[fc/n,(fc+l)/n], 1 
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(l[fc/n,(fc+l)/n], l[j/n,(j+l)/n])^(l[Vn,(/+l)/n], l[Vn,{i+l)/n])^(l[Vn,(/+l)/n]' 



n-1 



As in pLlj, when a, /?, 7, 5 are fixed, we can decompose the sum appearing in Br^a,i3,-y,sin) as 
fohows: 



E + E + E + E + E 

i=j=k=l i=j=k i=j=l i=k=l j=k=l 

\ i^i k^i ii^i i^j ) 



( 



+ 



E + E + E 



i=j,k=l i=k,j=l i=l,j=k 
\ k^i j^i j^i ) 



I 



+ 



+ 



E 



i=j,k^i i=k,j^i i=l,k^i j=k,k^i j=l,k^i k=l,kj^i 

\k^l,l^i m,k^l k^jjyti k^l,l^i k^l,l^i k^jj^i/ 



i,j,k,l 
are all distinct 



where the indices run over {0, 1 . . . , n— 1}. Similar computations as in [TT] show that the first, 
second and third sums are 0(l/(nlog^ n)); the fourth and fifth sums are 0(l/(nV9 log2 n)); 
the sixth, seventh and eihght sums are 0(l/(n^/'^ log^ n)); the ninth, tenth, eleventh, twelfth, 
thirteenth and fourteenth sums are also 0(l/(n^/'? log^ n)). We focus only on the last sum. 
Actually, it is precisely its contribution which will indicate the true order in ()1.6p . Once 
again, we split this sum into 24 sums: 



E + E 

k>l>i>j k>l>j>i 



(4.3) 



We first deal with the first one, for which we have 
1 



(n log n) 



^ PH{k - lY+'pH{i - JY'-'PHik - iTpHik - 3?ph{1 - iVPH{l - j)' 



k>l>i>j 



(n log 


n)2 


1 




(n log 


nY 


1 




\ 



k>l>i>j 

E E(^ - E(^ - 0-^'-^-^)/^ E(^ - ^r^^^^)/^ 

k l<k i<l j<i 



-{q-r-l)/q 



n log n 



n~l ^ n— 1 ^ 

'y^^-(r+l)/g ^ _^y^/-l^(r'+l)/g^l-(r+l)/9 _Ji 
^ ^nlog^n^ Incr 



logn 



where the notation a„ ^ 6„ means that sup„>i |an|/|6„| < +00. Since the other sums in 
(14. 3|] are similarly bounded, the fifteenth sum is 0(1/ logn). Consequently: 

9-2 

Y^E\Ar(nf] = 0(\l\ogn). 

r=0 

Finally, together with ()4.3p . we obtain £^[(1 - | ll-DS'nlll)^] = 0(l/logn) and the proof of 
(11. 6p is achieved thanks to Theorem 1 1.3[ 
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